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2.  Objectives 

Due  to  advaiict's  in  different  areas  of  microwave-  iuid  millimeter- wave  tedmiques, 
deimuid  for  efficient  CAD  tools  Inis  grown.  This  statement  is  today  ivs  true  iis  it  was  at 
our  U\st  workshop  held  two  years  ago  in  Munich  on  this  lojiic.  The  numerical  amvlysis 
in  time  donnun  is  attractive  since  it  describes  the  evolution  of  physical  quantities  in  a 
natural  way.  The  two  basic  concepts,  the  modelling  of  fields  and  networks,  .U'e  more 
closely  relaterl  than  it  seems  at  first  glance.  In  the  TLM  method,  for  examiile.  the 
Jietwork  model  is  the  basis  for  modelling  electromagnetic  fields.  Another  example  is 
the  application  of  genenilized  S-matrix  methods  and  diacojitir's  in  field  theory.  The 
purpose  of  this  amference  is  to  stimulate  synoptic  considerations  of  field  and  network 
theory  <uul  to  promote  a  lively  exduuige  betwwn  researchers  engagi'd  in  these  fields. 

3.  Venue 

The  workshop  will  take  |)lace  at  the  Hotel  Ambas.saflor  Berlin,  Bayreuther  Str.  42/43, 
10787  Berlin,  Germany,  Tel.:  +49/30/21902-0,  Fax:  +49/30/21902-380. 


4.  Conference  Secretary 

For  any  questions  or  infonnation  alxmt  the  workshop,  pleiuje  contact  the  Program 
Co-Chairpersons  or 

Michael  Krmnijholz,  Conference  Secretary 
Ferdinand-Braun-lnstitut  fiir  Hochstfrequenztechnik  Berlin 
Rudower  Chaussee  5,  12489  Berlin 
Tel:  +49/30/6392-2625,  Fax:  +49/30/6392-2612 
e-mail:  krumphol@fbh.wtza-berlin.de 


5.  Registration 

Advance  Registration  may  be  done  by  the  Registration  Form  enclosed  in  this  Prelimi¬ 
nary  Program.  On-site  registration  is  also  possible.  The  registration  office  at  the  Hotel 
Ambassador  will  be  open  at  October  28,  from  8:00  a.m.  to  10:30  p.m.. 

Please  note  that  late  registration  as  well  as  on-site  registration  is  subject  to  a  late  fee. 


6.  Projection  Facilities 

Projection  equipment  avixilable  will  be  a  standard  slide  and  an  overhead  projector. 

7.  Accommodation 

Blocks  of  rooms  have  been  reserved  at  the  following  hotels  near  or  at  the  workshop 
site  in  Berlin  at  special  workshop  rates.  Please  make  your  reservations  directly  with 
the  hotel.  Refer  to  the  Ferdimuid-Braun-lnstitut  fiir  Hochstfrequenztechnik  to  obtain 
the  special  rate. 

Hotel  Ambassador  Bayreuther  Str.  42/43 

Tel.:  +49/(0)30/21902-0  10787  Berlin 

Fax  :  +49/(0)30/21902-380  205  DM  single,  225  DM  double 

Hotel  Westerland  Knesebeckstr.  10 

Tel.:  +49/(0)30/3121004  10623  Berlin 

Fax  :  +49/(0)30/3136489  100  DM  single,  155  DM  double 

Hotel  Chariot  Giesebrechtstr.  17 

Tel.:  +49/(0)30/3234051  10629  Berlin 

Fax  :  +49/(0)30/3240819  60/90/115  DM  single,  170  DM  double 

Please  note:  The  rooms  are  only  reserved  until  September  14,  1993. 

8.  Transportation 

Berlin  is  served  by  three  airports:  Berlin-Tegel,  Berlin-Teini)elhof  ixnd  Berlin-Schone- 
feld.  Berlin-Tegel  is  situated  in  the  north-west  of  Berlin.  You  may  take  a  taxi  (about 
30  DM)  or  the  bus  no.  109  directly  to  Berlin  Zoologischer  G<\rten  (Berlin-Zoo)  (35 
min.).  Berlin-Temjjelhof  located  in  the  south-east  of  Berlin.  You  may  choose  a  taxi 
(alx)ut  25  DM)  or  the  underground  (U-Balm),  U8  to  Hallesches  Tor  and  change  to  the 
Ul  for  going  to  Berlin-Zoo.  If  you  arrive  at  Berlin-Sdionefeld  in  the  very  south-east  of 
Berlin,  the  most  jjreferable  way  to  get  into  the  center  is  taking  the  S-Balm  directly  to 
Berlin-Zoo  (70  min.).  A  taxi  will  be  about  60  DM.  In  general,  public  transportation  in 
Berlin  is  the  most  preferable  way  of  tnuisportation  within  the  city.  During  day-time, 
the  underground  and  buses  run  every  five  or  ten  minutes. 

If  you  are  cirriving  by  train,  the  trains  will  stop  at  Berlin-Zoo  or  Berlin-Hauptbahnhof 
from  where  you  can  take  the  S-Balm  directly  to  Berlin-Zoo.  If  you  are  arriving  by 
car,  leixve  the  Autobdmring  no.  10  at  Autolmhndreieck  Drewitz  to  the  Autobahn  no. 
115.  At  Autobahndreieck  Funktunn,  head  towards  north  (look  for  the  direction  to  the 
airport  Berlin-Tegel)  mid  get  off  the  motorway  at  Kaiserdamm  to  ixrrive  at  Berlin-Zoo 
via  Ernst-Reuter  Platz.  For  your  local  orientation,  we  have  enclosed  two  maps. 

9.  Program 

This  program  consists  of  the  invited  talks  and  the  non-iiivited  (oiitrilmtioii.''.  L.ii'- 
contributions  prcisented  in  .session  B‘2  will  he  accepted  at  the  wurk.sliop. 


Thursday,  October  28,  1993: 


8.00  Registration 


Session  Al 

8:30  Opening  Session 

8:40  Time  Domtun  Wave-Oriented  Data  Processing 

or:  How  to  Extract  Phenomenology  from  Observations 

9:25  Comparison  of  Different  Fie’d  Theoretical  Methods  of 
Ajialysis  of  Distributed  Microwave  Circuit  Elements 
10:10  -Break- 

10:30  Time  Donuun  Electromagnetic  Field  Computation 
with  Finite  Difference  Methods 
11:15  Finite  Difference  Time  Domain  Models  for  Coplamu' 
Waveguide  Discontinuities 

12:00  Enhanced  FDTD  Method  for  Active  tuid  Passive 
Microwave  Structures 
12:45  -Break- 


Session  A2 

13:45  On  the  Field  Theoretic  Foundation  of  the  TVimsmission 
Line  Matrix  Method 

14:30  TLM  Modelling  of  Guiding  and  Radiating  Microwave 
Structures 
15:15  -Break- 

15:30  Time  Domain  Simulation  of  Non-Linear  Networks 
Containing  Distributed  Interconnect  Structures 
16:15  Modelling  of  Planar  Microwave  Structures  in  Frequency 
and  Time  Domain 


Social  Event 

19:00  Piano  Recital:  Diana  Lawton 

Steinway-Haus,  Hardenbergstr.  12 
20:30  Berlin  Evening 

Berlin  Pavilion,  StraJSe  des  IT.Juni  100 


L.  Felsen 

R.  Sorrentino 

T.  W'eiland 

V.  Fouad  Hanna 

B.  Houslimand,  T.  Itoh 

P.  Russer,  M.  Krumpholz 

W. J.R.  Hoefer 

M. l.  Sobhy,  E.A.  Hosny 
R.  Valildieck 


iViday,  October  29,  1993: 


Session  B1 


8:00 

Dynamic  Simulation  of  Semiconductor  Devices 

R.  Steuzel,  W.  Klix 

8:45 

Signal-Processing  Approach  to  Robust  Time-Domain 
Modelling  of  Electromagnetic  Fields 

A.  Fettweis 

9:30 

-Break- 

9:50 

New  Results  in  Transient  Analysis  of  Crystal 

Ch.  Schmidt-Kreusel, 

Oszillators 

W.  Mathis 

10:10 

Adaptive  Detection  and  Tracking  of  Active  Scatterers 
by  Cascaded  Notch  Filters 

M.  Zouak,  J.  Saillard 

10:30 

Sub-mm  Wave  Circuit  Characterization  Using  the 
Finite  Difference  Time  Domain  Method 

N.I.  Dib,  L.P.B.  Katehi 

10:50 

FDTD  Modelling  of  Wirebond  Interconnects 

E.  Pillai,  C.  Bornkesse!, 
W.  Wiesbeck 

11:10 

Distributed  Computing  for  Transmission  Line  Matrix 
Method 

P.P.M.  So,  W.J.R.  Hoefer 

11:30 

TLM:  Order  of  Accuracy  Enhancement 

D.  de  Cogan,  A.  Soulos, 

P.  Enders 

11:50 

-Break- 

•  m 


Session  B2 


13:00  Space  and  Time  Discretization  in  Field  Computation 
Using  TLM 

13:20  Rigorous  and  Fast  Computation  of  Modal 
Johns’  Responses 

13:40  Transmission-Line  Matrix  Modelling  and  Huygens’ 
Principle  or  The  Range  of  Applicability  of  TLM 
14:00  Towards  Better  Understanding  of  the  SCN  TLM 
Method  for  Inhomogenous  Problems 
14:20  Generation  of  Lumped  Element  Models  of 
Distributed  Microwave  Circuits 
14:40  Towards  Exactly  Modelling  Open/ Absorbing 
Boundaries 
15:0C  -Break- 
15:20  Late  Contributions 

16:00  Concluding  Session:  open  forum,  panel  discussion, 
approximately  finished  by  17:00 


C.  Christopoulos, 

J.L.  Herring 

M.  Mongiardo,  M.  Righi, 

R.  Sorrentino,  W.J.R.  Hoefer 
P.Enders 


M.  Celuch-Marcysiak 


P.  Russer,  M.  Righi, 

C.  Eswarappa,  W.J.R.  Hoefer 
P.Enders,  A.J.  Wlodarczyk 
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Speciial 


Recital  Program 

Beriin,  Germany,  October  1993 


Diana  Lawton,  Piano 


I 

W.A.  Mozart  Sonate  d-dur,  KV.  311 

Allegro  con  spirito 
Andante  con  espressione 
Rondo 


C.  Debussy  L’Isle  Joyeuse 


Intermission 


F.  Chopin  Sonate  III  h-moll,  Op.  58 

Allegro  maestoso 
Scherzo  molto  vivace 
Largo 

Finale  presto,  ma  non  tanto 


1  Hotel  Ambassador 

2  Hotel  Westerland 

3  Hotel  Chariot 

4  Steinway-Haus 

5  Berlin  Pavilion 

6  Bahnhof  Zoo(-Iogisdier  Garten) 

7  Autobalmdreieck  Piinkturm 

8  Kaiserdainni 

9  Ernst-Reiiter  Platz 


Registration  Form 


Name/IEEE  number 


Address  . 

City  and  Zip-code  . 

Country  . 

Phone  ( . ) 

Fax  ( . ) 

1  wish  to  present  a 

late  contribution  on 

the  following  topic:  . 


Registration  Fee 

(includes  proceedings,  lunch  and  social  event) 

Please  note: 

Late  registration  as  well  eis  on-site  registration  is  subject  to  a  late  fee. 

IEEE  members  having  paid  by  1.10. 1993  $200/320DM  .1) . 

Non-members  having  paid  by  1.10.1993  $250/400DM  $ . 


IEEE  members  paying  aP'^r  1.10. 1993  $240/385DM  $. 

Non-members  paying  1.10.1993  $300/480DM  $, 

Extra  tickets  for  the  social  event 

number .  @  $45/75DM  $ 

total  enclosed  $. 


Please  make  checks  payable  to  the  Ferdinand-Braun-Institut  fiir  Hbchst' 
frequenztechnik  Berlin.  No  credit  cards  accepted. 


Remit  this  form  to:  Michael  Krumpholz,  Conference  Secretary 

Ferdinand-Braun-Institut  fur  Hochstfreciuenztechnik 
Rudower  Chaussee  5,  12489  Berlin,  Germany 


Time  Domain  Wave-Oriented  Data 
Processing 

OR 

How  TO  Extract  Phenomenology  from 
Observations 


L.  B.  Felsen 

Dept,  of  Electrical  Engineering 
Polytechnic  University  Six  Metrotech  C’enter 
Brooklyn,  NY  11201,  USA 


ContimuUly  increttsing  capabilities  for  sophisticated  measurement  and  nu- 
ineriatl  nuxlelling  are  now  available  to  generate  data  on  signal  transmission 
through,  and  scattering  by,  environments  of  subsUmtial  com])lexity.  Intelli¬ 
gent  extraction  of  information  from  the  received  signtd  for  surveillance,  iden¬ 
tification,  cU\ssification  and  other  tasks  is  aiderl  by  linking  features  in  data  to 
physical  features  in  the  enviroment  via  the  wave  processes  oiierative  under 
the  given  conditions.  To  implement  such  wave-oriented  data  i)rocessing,  it  is 
suggestive  to  assemble  a  dictionary  of  individual  scattering  scetiarios  to  de- 
tennine  the  footprints  established  by  each  in  the  observed  signal.  Data  from 
an  actuid  physical  or  numerictil  experiment  can  then  be  projected  onto  the 
dictionary  to  see  which  of  the  reference  signatures  are  conttiined  in  the  data 
base.  The  method  is  based  on  loc<d  forward  and  backward  tracking  of  wave- 
fields  in  the  phase  sijace  which  combines  both  configurational  (space-time) 
and  spectral  (wavenumber-fretiuency)  constituents.  The  data  i)roce.s.sing  is 
IJerfonned  via  windowed  transforms  and  multiresolution  teclmitiues.  Exmn- 
l)les  illustrate  frequency  mid  time  donuun  plu\se  space  tracks  jiroduced  by  a 
variety  of  basic  scattering  sceiuu-ios,  with  numericid  examijles. 


Signal- Processing  Approach  to 
Robust  Time-Domain  Modelling  of 
Electromagnetic  Fields 


Alfred  Fettweis 

Ruhr-Universitat  Bochum 
Lehrstuhl  fiir  Nachrichtentechnik 
44780  Bochum,  Germany 


In  recent  publications  jl]  -  (6|,  a  new  method  for  integrating  piU'tial  dif¬ 
ferential  e<|uations  describing  ivctiud  physicjU  systems  has  been  presented. 
This  metho<l  is  based  on  simulating  the  iictual  continuo\is-dom<un  system 
by  means  of  a  discrete-domain  system,  and  this  in  such  a  way  that  the  fol¬ 
lowing  features  hold: 

1.  Preservation  of  origimilly  existing  passivity  and  incremental  passivity, 
and  this  in  such  a  way  that  these  pro|)erties  become  available  in  the 
multidimensional  (MD)  scMise  even  though  they  existed  originally  only 
in  the  one-dimension<\l  (1-D)  stnuse  (i.e.  with  rv-.spect  to  time).  As 
a  result,  one  can  <vchieve  not  only  full  stability  with  respect  to  the 
discretiziUion  in  si)<ice  and  time  but  also  full  stability,  and,  more  gen¬ 
erally,  full  robustness  with  respect  to  the  computational  errors  that  are 
due  to  rounding/tnmcation  iuul  overflow  corrections  and  to  extrcuieous 
sources. 

2.  Preservation  of  the  massiv  ])andlelism  and  the  exclusively  local  nature 
of  the  interconnections,  which  is  inherent  to  all  physical  systems  with 
finite  propagation  speed.  As  a  result,  for  any  given  fixed  time  instant  to 
l)e  considered,  the  computations  c<m  be  Cixrried  out  .simultaneously,  thus 
fully  in  paniil 'I,  m  all  the  spatial  .s<unpling  points,  and  the  computations 
in  any  of  these  points  reciuire  previously  computed  results  only  from 
the  immefliate  neighbouring  points. 

3.  Arbitrarily  changing  i)arameters  as  well  as  a  arbitrary  buundary  sha])es 
and  conditions  can  b«?  taken  into  <\ccount  in  a  .straightforward  manner. 

In  order  to  lU'hieve  recursibility  (amiputability),  the  simulation  may  not  lx- 
based  on  the  field  variables  ajiiiearing  in  the  original  partial  different  ial  eciwa- 
tions.  Instead,  mrnsponding  .so  called  wav»;  variables  shouhl  lie  emiiloyed. 


thus  vtiriables  of  the  type  occuriug  in  relation  with  the  scattering-matrix  for¬ 
malism.  This  way,  the  mecluuiism  involvetl  in  the  physical  system  becomes 
inteqjretable  tui  incidence-to-scattering  (reflection,  tnmsmissior.)  mecha¬ 
nism,  i.e.  a  mechtuiism  exhibiting  a  cause-to-eflect  (causality)  relationship. 
The  latter  in  turn  gives  rise  to  computational  rules  that  exhibit  the  sequential 
nature  needed  for  obtaining  an  algorithm. 

It  appetu-s  easiest  to  ajiidy  the  method  by  first  representing  the  system  by 
means  of  a  multidimensional  electric  circuit.  From  this,  the  desired  algorithm 
ciui  be  derived  i)y  aiJjjlying  the  staiKhu-d  procedures  know  from  the  theory  of 
multidimensional  wave  digital  filters  (7|,  which  has  originally  been  developed 
within  the  context  of  digital  signal  processing.  The  ai)proiU'h  is  api)licable 
without  difficulty  to  systems  described  by  Maxwell’s  etjuations. 

References 
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Finite  Difference  Time  Domain  Models 

FOR  COPLANAR  WAVEGUIDE  DISCONTINUITIES 


Victor  Found  Hanna 


FVance  Telecom,  Centre  National  d’Ctudes  des  Telecommunications 
Centre  Paris-B,  38-40  rue  du  General  Leclerc 
92131  ISSY-LES-MOULINEAUX,  FRANCE 


This  paper  assembles  various  moclels  that  were  developped  using  the  Finite 
Difference  Time  Domain  (FDTD)  Method  to  characterise  cojjlanar  waveguide 
(CPW)  uniaxial  discontinuities  (step  in  widths  gap,  short-circuit,  open  cir¬ 
cuit)  and  multimaxial  ones  (T  junction,  bends).  Comparison  of  these  models 
with  those  develoi)ped  using  other  general  puqwse  time  tlomain  methods  like 
the  TLM  method  or  frequency  dommn  ones  like  the  finite  element  method 
will  be  presented  whenever  i)ossible.  Models  for  other  possible  CPW  disconti¬ 
nuities  use<l  for  jnonolithic  api)lications  like  mr  bridges,  connections  between 
lines  on  different  substrates  and  via  holes  for  backed  grounded  CPW,  will  be 
given  also.  For  CAD  purposes,  each  of  these  discontinuities  will  be  modelled 
between  predefined  reference  i)Uuies  either  by  an  etiuivalent  electric  circuit 
or  a  scattering  matrix. 


TLM  Modelling  of  Guiding  and  Radiating 
Microwave  Structures 


Wolfgaiig  J.  Ft.  Hoefer 

NSERC/MPR  Teltech  Research  Ch.\ir  in  RF  Engineering 
Dept,  of  Electrical  and  Computer  Engineering,  University  of  Victoria 
Victoria,  B.C.,  Canada  V8W  3P6 


Tliis  p»per  provides  an  overview  of  recent  progress  in  modelling  and  simula¬ 
tion  Oi  v-iectromagnetic  structures  with  the  TLM  method.  It  is  divided  into 
three  pjtrts:  New  theoretical  developments,  computer  imiilementation  and 
validation,  imd  design/optimisation. 

New  theoretical  develojunents  include  generalized  rectangular  mesh  al¬ 
gorithms  for  increased  flexibility  of  discretization,  new  multi-modal  Johns 
matrix  techniques  for  high-quality  absorbing  boundary  modelling  in  sei)a- 
rable  structures,  and  higher-order  one-way  equation  models  of  absorbing 
boundaries  for  dispersive  guiding  structures  as  well  as  for  open  radiating 
structures,  such  as  microstriiJ  patch  iuitennas.  Methods  for  second-order 
accurate  di.scretization  of  non-cartesian  geometries  are  .dso  described  in  the 
first  part. 

Computer  implementation  and  validation  of  TLM  algorithms  are  impor¬ 
tant  asiiects  that  detennine  the  particaJ  usefulness  and  reliability  of  the  TLM 
method.  The  second  part  of  this  i)aper  thus  comjjares  the  serial  and  par¬ 
allel  progrmnm  structures  of  2D  and  3D  TLM  simulators,  describes  their 
performance  on  tyijical  comi)uter  platforms,  and  conveys  a  flavour  of  the 
simulator-user  interface.  Solutions  of  sevenJ  canonictd  jiroblems  demon¬ 
strate  the  accunicy  of  TLM  results  for  time  and  frec|uency  domain  charcu;te- 
ristics  of  microwave  components,  as  well  as  for  visualizing  fields  and  circuit 
functions. 

Design  and  optimisation  are  the  ultimate  engineering  ap])lications  of  a 
modelling  technique.  Two  possible  api>roaches  are  described  in  the  third 
part.  The  first,  more  chissical  approtich  combines  time  domain  field  analysis 
with  jui  optimization  program.  By  using  datapi])e  teclmicjueH,  the  TLM 
analysis  can  be  perfonned  on  a  powerful  i)roces.sor  siich  as  .MAS PAR  or  the 
Coimection  Machine,  or  the  task  can  be  divided  cunong  several  workstat bjn.s. 
Tile  other  approach  makes  iLse  of  time  reversiil  to  syntht^size  a  structure 
geometry  directly  from  its  desired  frequency  response.  This  exciting  new 
procedure  is  explained  and  demon.straterl  irsing  a  .specific  syntlnwis  exaini)le. 

The  paper  concludes  with  a  look  at  |M)s.siljilities  and  <level()pinents. 
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It  is  well  known  that  the  FDTD  (Finite  Difference  Time  Domain)  method 
is  a  powerful  method  for  analyzing  the  electromagnetic  wave  behaviors  in  a 
complicated  geometry.  However,  it  is  a  memory  ijitensive  and  time  consu¬ 
ming  o|>eration.  Recently,  the  author’s  group  has  invented  and  imjjlemented 
several  techniques  to  alleviate  these  deficiencies.  Sijecifiailly,  we  have  im¬ 
plemented  the  FDTD  Diakoptics  methods  to  use  numerical  Green’s  function 
to  replace  large  comiJUtation  volume  with  its  impulse  resjjonse.  Hence,  the 
memory  re<iuirement  is  drastically  re<luced.  For  ivcceleration  of  numerical 
speed,  we  have  implemented  a  method  ba^d  on  the  system  identification 
technique.  A  reduction  computation  time  of  a  factor  of  ten  cfui  readily  been 
attained.  In  addition,  for  the  first  time  in  the  electromagnetics  community  in 
the  world,  we  c<ui  deal  with  a  large  volume  containing  active  and  nonlinear 
regions  (or  devices,  phisma,  etc.)  by  means  of  FDTD  environment. 

In  this  talk,  fundamental  conceiJts  and  efforts  for  their  implementations 
are  presented.  Several  examples  are  disaissed.  Since  the  methorls  are  still 
in  the  development  stage,  some  of  the  recent  findings  are  also  added  at  the 
time  of  presentation. 
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The  TLM  method,  developetl  and  first  published  in  1971  liy  Johns  tuid  Beurle 
has  emerged  as  a  ix)werful  method  for  computer  modelling  of  electromagnetic 
fields.  In  TLM  the  space  is  subdivided  into  cells.  The  electromagnetic  field 
dynamics  is  modelled  by  wtwe  pulses  proi)agating  between  adjacent  cells  and 
being  scattered  within  the  cells.  The  main  advantage  of  the  TLM  simula¬ 
tion  resides  in  the  capability  to  model  circuits  of  arbitrary  geometry,  and  to 
compute  and  to  display  the  time  evolution  of  the  fields.  The  TLM  method 
exhibits  an  excellent  numericiil  sUibility  and  is  also  suitable  for  modelling  of 
lossy,  disjjersive  iuul  nonlinear  me<lia. 

Introducing  the  Hilbert  sj^we  representation  for  the  field  state,  the  de¬ 
scription  of  geometrical  structures  and  the  field  evolution  is  ]7erformed  alge¬ 
braically.  Field  theoretical  foundations  of  the  two-dimensional  TLM  method 
and  the  three-dimensional  TLM  method  with  condensed  node  are  given  using 
the  Method  of  Moments  with  sectional  base  functions.  It  is  shown  that  the 
sampling  of  the  tangential  electric  mid  magnetic  field  conqionents  in  the  cell 
boundary  surfaces  yields  a  correct  bijective  mapjjing  between  electronif^ne- 
tic  field  components  and  TLM  wave  amplitudes. 

Tlie  space  discretization  with  reguhu-  meshs  as  well  with  nonuniform  and 
cvu*ved  meshs  and  the  error  introduced  by  the  mesh  discretization  me  discus¬ 
sed.  For  the  calailation  of  the  TLM  and  FDTD  disj)ersion  relations,  we  luse 
a  new  generalizetl  method.  A  criticid  comiJarison  of  two-dimensional  and 
three-dimensioiud  TLM  with  other  finite  difference  time  domain  methods  is 
given. 


Time  Domain  Simulation  of  Non-Linear 
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lYaiisniission  proi)erties  of  interconnects  such  as  signal  delay,  reflection,  at¬ 
tenuation,  dispersion  and  crosstalk  must  lie  taken  into  consideration  in  the 
analysis  and  design  procedures  of  high  speed  digitcU  circuits  and  high  fre¬ 
quency  microwave  systems,  it  is  important  to  include  the  transmission  line 
behaviour  of  interconnects  between  the  system  comijonents,  if  the  behaviour 
of  the  overall  system  is  to  lie  accurately  simulated. 

In  general,  the  system  can  be  divided  into  lumi>ed,  non-linear  subnet¬ 
works  and  tUstributed  interconnect  structures.  The  distributed  interconnects 
can  be  represented  by  different  biu?ic  models  of  transmission  lines,  uncoup¬ 
led  lossless,  coupled  lossless,  luicoupled  lossy  and  cou])led  lossy  trmismission 
lines. 

In  this  i)aper,  different  approaches  are  given  to  characterize  the  inter¬ 
connects  using  time  domain  or  frequency  domain  jjrocedures.  Each  inter¬ 
connect  structure  can  be  representetl  by  a  circuit  block  described  by  its  ter¬ 
minal  relations.  In  the  simple  case  if  the  interconnect  is  modelled  as  a  lossless 
and  non-dispersive  transmission  line  structure,  the  termiiud  relations  can  lie 
represented  by  a  set  of  difference  e<|uations.  In  id!  other  cases  the  terminal 
relations  can  l>e  derived  by  using  CiUculaterl  or  me<isured  .scattering  para¬ 
meters  of  the  interconnect  stntcture  in  the  time  domain  .s(f.).  which  can  be 
generated  directly  in  the  time  domain  or  by  transforming  frecpiency  domain 
data  to  the  time  donuun. 

If  the  interconnect  is  modelled  as  a  lossy  ciuasistatic  transmission  line, 
nodal  time  domain  analysis  can  be  applied  to  obtmn  directly  .s(<).  However, 
in  the  Civse  of  freciuency  dependent  transmi.ssion  line  parameters  and  when 


coupling  between  the  transinissiun  lines  is  not  litnitecl  to  acljiicent  lines,  the 
time  domain  nodal  approacli  does  not  apply.  In  the  case  of  lossy  and  disper¬ 
sive  transmission  line  structures,  s(t)  can  Ije  obtained  directly  by  using  time 
domain  full  wave  anylysis  (TLM  and  FDTD)  or  frequency  domain  analysis 
to  obtain  frequency  domain  scattering  iwiuneters  which  can  be  transformed 
to  obtain  s{t). 

Finally,  we  have  now  a  set  of  luinjjed,  non-linetvr  subnetworks,  all  are 
modelled  by  their  time  dontain  tenninal  relations.  A  general  procedure  based 
on  the  state  space  aiiproadi  has  been  developed  to  obtain  the  time  domain 
analysis  of  whole  netvrark. 

The  transient  responses  of  many  examples  have  l>een  obtiiined  by  using 
the  proposed  analysis  procedure.  These  examples  include  tUfferent  types 
of  intercoimect  structures  in  high  s])ee<l  digitd  circuits  and  high  frequency 
analogue  circuits.  The  advantages  of  the  proposed  methorl  <xre  the  high 
computational  stability  and  efficient  numerical  accuracy. 
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Numerical  methods  for  the  analysis  of  microwave  structures  belong  basically 
to  two  categories  corresponding  to  the  numerioil  formulation  of  Maxwell’s 
equations  in  differential  or  integral  form.  Strictly  munerical  methods,  on  the 
one  hand,  are  based  on  the  discretization  of  Maxwell’s  ecjuations,  (e.g.  Finite 
Clement,  Finite  Difference  (in  frequency  or  time  domain))  or  on  the  imple¬ 
mentation  of  Huygen’s  principle  (TVansmission  Line  Matrix  (TLM)  method). 
Tluinks  to  their  flexibility,  they  are  well  suited  for  i)roblems  with  irregular 
geometries.  On  the  other  hand,  methods  requiring  a  degree  of  mathematical 
preprocessing,  such  as  Spectnd  Domain  Approach  (SDA)  and  mode  matching 
(MM),  lead  to  very  efficient  computer  codes,  but  are  limited  to  problems  with 
simple  geometries. 

To  illustrate  the  aljove  asi)ects,  a  comparison  between  the  Finite  Dif¬ 
ference  Time  Domain  and  the  Mode  Matching  methods  as  ai}i)lied  to  the 
analysis  of  ))ractial  microwave  structures  will  Ije  presented  in  this  jjaper. 
The  methods  tire  compared  in  tenns  of  CPU,  memory  recjuirements,  versati¬ 
lity,  accuracy  of  the  results,  etc.  A  number  of  examples  of  components  and 
discontinuities  both  in  conventional  as  well  as  integrated  technologies  will  be 
discussed. 
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DeveloiJinent  and  optimization  of  novel  semiconductor  devices  demand  the 
use  of  two  (2D)-  and  three  (3D)-dimensional  simulation  methods. 

At  first  a  survey  of  the  different  jnodel  levels  of  device  simulation  (micros¬ 
copic  and  macroscopic  models)  will  be  given.  Subsetiuently  the  2D-simulator 
SEMICO  and  the  3D-simulator  SIMBA  developed  at  the  Dresden  Polytech¬ 
nic  University  of  Teclmoloy  and  Economics  juid  the  Dresden  University  of 
Technology  will  be  jiresented.  In  these  programs  drifl-diffusion  models  are 
used  and  the  semiconductor  e(|uations  are  solvetl  with  finite  difference  me¬ 
thods. 

Tlie  pajier  contiiins  a  discussion  of: 

•  basic  semiconductor  equations 

•  models  for  mobilities,  recombination,  generation,  temiierature  dej)en- 
denoe 

•  material  parameters 

•  boundary  conditions 

•  munerical  methods 

•  si>atial  discretization 

•  time  integration 

•  treatment  of  large  linear  systems  of  equations 

•  calculation  of  Y-i>arameters 

•  calculation  of  typical  jKiwer  giuns  and  of  sinall-  signal  e<iuivalent.  circuit 
motlels 

Tlie  different  simulation  possibilities  are  demonstrated  in  connection  wit  li 
novel  lliyV-heteroiuuciion  devices  fUBT-  HF.MT  IPtn  _ 
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Time-  and  frequency-domain  technicjues  are  complementary  tools  in  the  euia- 
lysis  and  design  of  microwave  and  oi)tical  integrated  circuits.  While  time- 
domain  metliods  are  very  aseful  for  wideband  itpplications  like  transient  i)ro- 
blems  as  well  as  for  nonlinejir  iuntlysis,  frequency-domain  techniques  tve  com¬ 
putationally  more  efficient  for  the  rnajority  of  microwave  design  jiroblems, 
which  are  of  comparatively  narrow  band  nature. 

To  address  the  different  design  jiroblems  within  the  framework  of  one 
analysis  techniqvie,  the  finite  difference  (FD)  technicine  in  the  frecpiency- 
domain  and  the  finite  difference  tiine-<lomain  (TDFD)  technicpie  c\re  availa¬ 
ble  as  complementary  tools.  For  the  tiine-tlonnnn  triuisiuissioii  line  matrix 
(TDTLM)  method,  no  freciuency-doinain  counter])art  existed  until  recently 
when  the  TLM  node  was  tested  in  the  frec|uency-domiun  in  conjunction  with 
a  novel  S-parameter  extraction  technic|ue.  The  resulting  frecjuency-tlomiun 
TLM  (FDTLM)  technique  is  a  very  pr)werful  and  flexible  numerical  modelling 
tool  for  frequency-domain  design  jiroblems.  Its  computational  efficiency  and 
flexibility  is  better  than  most  fre<|uency-domajn  technic|ues  known  today. 
Since  the  siiace  discretiziition  is  using  the  same  transmission  line  nodes  as  in 
the  TDTLM,  the  FDTLM  re])resents  a  true  fre<iuency-<lomain  counteri)i\rt 
to  the  time-domain  TLM. 

In  this  paper  a  general  overview  will  lie  given  about  modelling  microwave 
ciraiits  in  the  time-  and  frequency-domaiii.  New  research  results  concerning 
the  FDTLM  will  lie  presented  mid  its  relationshij)  to  the  TDTLM  will  be 
discussed. 


domain.  This  is  a  ty])ical  example  wliere  l)oth  fre(|ueucy  tuid  time  domain 
analysis  are  essential  and  only  the  combination  yields  the  successful  result. 

Typical  engineers  may  wonder  why  one  applies  time  domain  analysis  to 
basically  monochromatic  held  problems  at  lUl.  The  anwser  is  simple:  it 
is  mudi  faster,  needs  less  computer  memory,  is  more  general  and  typically 
more  accurate.  Sjjeed  up  hictors  of  over  200  have  been  reached  for  retUistic 
problems  in  filter  and  wave  guide  design. 

The  small  core  s|)ace  requir»neiit  makes  time  domain  methods  a)}plicable 
on  desktop  computers  using  millions  of  cells,  and  six  unknowns  per  cell  -  a 
dimension  that  has  not  yet  been  reached  by  frequency  domiiin  appixuvdies. 
This  enormous  amount  of  mesh  cells  is  absolutely  necessi\ry  when  complex 
structures  or  structures  with  spatial  dimensions  of  many  wavelengtlis  are  to 
be  studied.  Our  personal  record  so  far  is  ti  waveguide  ])rublem  in  which  we 
used  72.000.000  miknowns. 

We  will  present  a  number  of  conventional  examples  from  electrical  en¬ 
gineering  and  the  forefront  research  in  high  ixjwer  RF  generation  to  demon¬ 
strate  that  time  domain  methods  are  ’Ijetter  by  design’  and  that  there  is  no 
pure  frequency  domain  method  around  tluU  ciui  reach  the  smue  generality 
and  practicability. 
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The  solution  of  Mivcwell’s  eciuations  in  time  domain  has  now  been  used  for 
more  thiui  two  deawles  and  has  had  great  success  in  imuiy  different  appli¬ 
cations.  The  main  attraction  of  the  time  domain  approach,  originating  in  a 
paper  of  Yee,  is  its  simplicity.  It  takes  only  marginal  effort  to  write  a  compu¬ 
ter  code  for  solving  a  simi)le  scattering  problem  compared  with  conventioiuil 
frequency  domain  methods. 

However,  when  applying  the  time  domain  ajjproach  in  a  gimeral  way  to 
arbitrarily  complex  jiroblems,  many  se<;mingly  simjjle  additioiiid  ])roblems 
arid  up. 

We  present  here  a  theoretical  framework  for  solving  Maxwell’s  equations 
in  integral  from,  re.sulting  in  a  set  of  matrix  ecpiations,  each  of  which  is  the 
discrete  analogon  to  one  of  the  original  Miixwell  etpiatitnis.  This  api^roach  is 
called  Finite  integration  Theory  and  was  first  deveIoi)ed  for  frec|uency  domain 
problems  also  about  two  decades  i\go.  The  key  point  in  this  formulation  is 
that  it  can  lie  a))plied  to  static,  harmoiiic  and  time  dej)endent  fiekls,  mainly 
because  it  is  nothing  but  a  computer-compatible  re-fonnulation  of  Mcxxwell’s 
equations  in  integral  form. 

When  s])ecii\lized  to  time  domain  fields,  the  method  actually  contains 
Yee's  algorithm  as  a  sub.set.  Further  iulditions  include  los.sy  materials  and 
fields  of  moving  charges,  even  including  fully  relativistic  analysis. 

For  many  pn\ctical  problems  the  pure  time  domiun  algorithm  is  not  suf¬ 
ficient.  For  instance  a  waveguide  transition  analysis  rec|iiires  the  knowledge 
of  the  incoming  and  outgoing  nuKle  i>atterns  for  i)roi)er  excitation  in  time 


